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Abstract 

A master equation for the evolution of two-dimensional universe 
is derived based on the simplicial quantum gravity regarding the evo- 
lution as the Markov process of a space-time lattice. Three typical 
phases, expanding, elongating and collapsing phase, which have been 
found in the numerical simulation, are studied together with their 
boundaries, analytically Asymptotic solutions of the evolution equa- 
tion for statistical quantities, such as the averaged area, the boundary 
length, correlations of their fluctuations, are obtained for each phase 
and boundary. After introducing a physical time the cosmological 
significance of each phase is discussed. 

PACS: 98.80.Qc, 04.60.-m, 98.80.Bp, 02.50.Ga 



1 Introduction 



Recent observations and their analyses on anisotropies of the cosmic mi- 
crowave background [IJ [2] have revealed the large scale correlation between 
two points separated far beyond the event horizon. It is considered as a relic 
of quantum fluctuation produced during the period when the universe is a 
quantum mechanical size. The full sky map of universe tempts us to imagine 
that we are glancing at the stage where the universe begins to explode. In 
order to describe this period on the scientific ground it is inevitable to have 
a tool which is capable of handling the creation of space, time, and matter 
all at the same time. 

It is quantum gravity which is aimed to describe such a circumstance. 
Quantum gravity is also expected to be the theory for quantizing the gravi- 
tational field. However, incompleteness of various theoretical attempts in the 
4-dimensional space-time by standard quantization methods made most of 
the investigations turn toward a totally new direction, i.e. the super-string 
theory, except a few investigations [31 H] . On the other hand, as far as the cre- 
ation of space, time, and matter is concerned, there exist consistent theories 
of quantum gravity in the 2-dimensional space-time, where no quantization 
of the gravitational field is necessary. It is our aim of this article to give 
an insight on birth and growth of the universe through the 2-dimensional 
quantum gravity. 

There are several standard methods of the 2-dimensional quantum grav- 
ity, for examples, the Liouville theoryj5] based on the standard field the- 
ory, the matrix model[6] representing the combinatorial geometry [7], and the 
dynamical triangulation[8j as its numerical realization. Recently, we have 
proposed a scenario of the initial stage of universe in terms of the simplicial 
quantum gravity [9j [TP] . With a few fundamental axioms a d-dimensional 
universe has been constructed as a set of complexes of d-simplices. Numer- 
ical simulations have been performed in 2, and 4 dimensions by applying 
techniques of the dynamical triangulation. They have indicated that there 
are basically three types of universes, which have been called the collapsing, 
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the expanding and the elongating phase. For the 2- dimensional case nu- 
merical results have been compared to predictions of the boundary Liouville 
theorypj], and shown to reproduce the expected results. 

However, those comparisons have been limited within the collapsing phase, 
where the stable vacuum exists and the standard field theory technique is 
applicable. In order to analyze non-static phases, in which the spaces keep 
growing indefinitely, we need to find an appropriate method that can handle 
a system without the stable vacuum. In this article we propose a method to 
treat such cases based on techniques of the non-equilibrium statistical me- 
chanics. This study shows that the 2-dimensional simplicial quantum gravity 
enables to describe most of the essential features of the 4-dimensional uni- 
verse, such as the inflation and the big bang, even under the limited condition 
of dimensionality. 

There exist a few attempts for creation of universe in the similar stand- 
point as ours, expecting the underlying discrete nature of space-time, such as 
the causal dynamical triangulation (CDT)[12J and the loop quantum gravity 
(LQG)[13]. When one consider the very beginning of universe, one inevitably 
give the initial hypothesis which cannot be proven. The main difference of 
those models comes from a point where the model starts. Both the CDT and 
the LQG assume the existence of the Lorentzian space-time and quantum 
dynamics before the birth of universe, while we consider all the physical laws 
appear at or after the birth of universe. Thus our future task will be to derive 
the dynamical law and the Lorentzian invariance of dynamics. 

In the next section (Sec. 2) the scenario of creation of universe is reviewed 
within the 2-dimensional quantum gravity. I derive a master equation as 
the continuous limit of the dynamical triangulation with a boundary. In 
Sec. 3, I study the statistical property of static universe through calculat- 
ing the partition function as a generating function for the average of area, 
boundary length, and correlation functions of fluctuations. Sec. 4 is devoted 
to overview possible non-static solutions of the master equation. In Sec. 5, I 
solve equations of the motion for averaged quantities derived from the mas- 
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ter equation for the expanding phase and its boundaries. Following section 
(Sec. 6) describes the elongating phase, where the boundary length stays con- 
stant while the area keeps increasing. In Sec. 7 two special states of universe 
are discussed. One is a universe with the minimum boundary length, which 
may be called the polymer state, appearing for a specific value of the param- 
eter related to the matter central charge, and the other is a boundary state 
along the line locating between the expanding and elongating phases. Sec. 8 
is devoted to discuss evolution of the universe in 1 + 1 dimension by intro- 
ducing a physical time, and show that it is possible to interpret each phase 
as a state of actual universe at various stages of its evolution, such as the 
inflation, the big bang, the stationary, and the crunch. Last section (Sec. 9) 
summarizes my aim, accomplishments, and basic problems left to be solved. 
Three appendices are in order for the master equation in the continuous limit 
and two basic equations of the motion referred in the text. 



2 Master equation for the evolution of uni- 
verse 

We define the 2-dimensional universe to be a set of all the possible simpli- 
cial complexes of triangles with one boundary. We shall call a complex as a 
configuration, or sometimes a world. It is possible to classify those simplicial 
manifolds by the number of handles g. In this article we limit our worlds 
only to those manifolds of the simple disk topology without handles, for sim- 
plicity. Each complex has its neighbors which can be obtained by attaching 
or detaching one unit triangle. In 2-dimension, there are, in general, four 
directions of neighbors for each configuration, either increasing or decreas- 
ing the area and boundary length by one unit. These procedures are called 
moves. 

In practice, a series of configurations is generated by selecting randomly 
one of those possible neighbors connected by a move at each step. We start 
with any but normally the simplest configuration made of one elementary 
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triangle to create a Markov chain of the length N, 

dii — > a 2 — > ■ ■ ■ — > ajy. 

A configuration is specified by a set of 'quantum' numbers, a = (N 2 , N\, 
k), namely, the total number of elementary triangles, N 2 , the number of 
edges on the boundary, N ly and the quantum number k identifying a specific 
diagram with a fixed (N 2 , Ni). A selected move is accepted if the resultant 
configuration satisfies following manifold conditions: 1) at most two triangles 
attach through one edge, and 2) triangles sharing one vertex form a 2-disk, 
or a semi-disk if the vertex is on the boundary 

Mathematically the dynamical triangulation can be described by a master 
equation for the probability p a [n] of a configuration a at the n-th step, which 
is expressed as 

p a [n + 1] = p a [n] + ^{pp[n]wp^, a -p a [n]w a ^}, (2.1) 

0(a) 

where the sum over (3(a) is understood to take all those neighbors {/?} con- 
nected by one of the four moves from a configuration a . The transition 
probability is chosen so that in the case when the equilibrium distri- 

bution p a [equi] exists, it satisfies the equation, 

p a [equi)w a ->/3 = pp[equi]wp^ a , (2.2) 

known as the detailed balance. From one of the axioms, which states that all 
the possible worlds are equally probable 0, Pa[eq ui] should be equal for any 
distinct configuration a. 

The equal probability axiom results the state independence of the transi- 
tion probability w a -+p between neighboring two configurations, and the mas- 
ter equation can be rewritten within the reduced space of the first two quan- 
tum numbers, a = (N 2 , Ni), as 

p a [n + 1] = p a [n] + ^{pblnlwb^a - p a [n]w a ^ b }, (2.3) 

b 

1 We shall call it the equal probability axiom. 
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where the sum over b runs from 1 to 4 corresponding to the 4 moves, which 
shift {N 2 ,Ni} to {iV 2 ± 1,N\ ± 1}. The detailed balance condition in the 
reduced space is then written as, 

pa[equi]w a ^ b = p b [equi]w b ^ a . (2.4) 

As a consequence of the equal probability axiom, p a [equi] should be propor- 
tional to the number of configurations, <ft a . The detailed balance condition 
can be satisfied by taking the transition probability w a ->b as 

<Pb , _v 

<Pa + <Pb 

It may be possible to choose another form for u) a _>6 which fulfill the detailed 
balance condition. Our specific choice of w a -^b will be justified when we take 
the continuous limit of the master equation. 

The master equation in the continuum space is obtained by taking the 
limit that the scale parameters {A , / , To}, introduced for the area A = A N 2 , 
the boundary length, I = IqNi, and the Markov process tinted, r = Ton, tend 
to zero. Leaving details of the derivation in Appendix A, it is written as 

d r , , ! r / d 2 d dS a \ r d 2 d dS a \^ r . , ■, 

^ w = { W + MJX ) + (a? + « If) M A < T ^ < 2 ' 6 > 

where S a = 5(^4, /) is defined by 4> a = exp(— S a ). Scale parameters are tuned 
as 

A 2 I 2 

Eq.(2.6) clearly shows that the distribution a is the stationary solution, 
and it is the reason of our choice for the transition probability, w a ->b- The 
function cf) a was computed numerically by the dynamical triangulation with a 
boundary |14j. and later the exponent was analytically given by the boundary 
Loiuville theory pTj and the matrix theory [TO] as 

I 2 

S(A, I) — -alogA — (3 log l + j-+fj,A + /A (2.7) 



2 We may call it simply time unless we need to distinguish it from the physical time 
which will be defined later. 
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where a = —1 — 1/b 2 , j3 = —2 + 1/b 2 , and 7 = l/(4sin7r& 2 ). The parameter 
b is related to the Liouville background charge Q through Q = b + 1/b. In 
order to avoid the singular behavior of the master equation at those b where 
7 diverges, we rescale the area, the boundary length, and the time by 7, as 
A = 2/7, I = y/7, and r = t/7 2 , respectively, which leads 

-P\x, y;t) = i— + ± dS ^ y ^ +^ + ^- dS ^ y) \ p \ x , y t], (2.8) 
dt •> dx 2 dx dx dy 2 dy dy > 

with 

v 2 

S(x,y) = -cdogx - logy + — + \x + \ B y + 31og7. (2.9) 

x 

The cosmological constant, A, and the boundary cosmological constant, X B , 
are also rescaled accordingly as /i = 7 A and \i B = jX 8 . The last term in 
S(x, y) modifies the overall normalization of the distribution, and we shall 
put this term in the normalization constant. 

In the two dimensional parameter space (A, X B ) the numerical simulation 
for the pure gravity {1/b 2 = 3/2) has revealed three critical lines separating 
into regions corresponding to three distinguished phases of universe. We 
shall write each phase as PI, P2, and P3, and name them as the collapsing, 
the expanding, and the elongating phase, respectively. We also write the 
boundaries between PI and P2 as Bl, P2 and P3 as B2, and PI and P3 as 
B3. The relative locations of phases and boundaries are depicted in Figure^] 
The numerical simulation has shown that the area tends to decrease towards a 
point universe when the parameters (A, X B ) are in the region of the collapsing 
phase, which locates most of the right half-plane of the parameter space, while 
the area grows indefinitely with the maximum rate when the parameters are 
outside of PI. This non-static region is divided into two phases by the line 
X B = in the A < half-plain. Above the critical line the boundary length 
tends to a constant either like an elongating tube or string with the minimum 
radius. Below the critical line both the area and the boundary length increase 
with the maximum rates. 

In the subsequent sections we shall obtain possible solutions for various 
phases, which are classified according to the behavior of average area and 
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Figure 1: Collapsing phase(Pi). Expanding phase(P2). Elongating phase(P5). 
The boundary state is on the boundary B2, and the Polymer state is in the elon- 
gating phase P3. 

boundary length in the parameter space (A, X B , /3). 

3 The collapsing phase (PI ) 

We begin with studying the static solution p a [equi] ~ exp(— S a ) of the master 
equation. Integrals in the generating function, 




(3.1) 
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are in general diverging, and we need to regularize it. We regard it as the 
continuum limit of the matrix model partition function, 

Z M [^fi B }= £ exp(-/iiV 2 - // B iVi) Z M [N 2 ,N X ], (3.2) 

(N 2 ,Ni) 

where the sums over (N 2 , Ni) are restricted by N 2 > N% — 2, N\ > 3, and 
N2 + N\ = even integer, due to the manifold conditions we have employed. 

The number of distinct congigurations has been calculated either by the 
graph theory [7J or the matrix mo del [B], which approaches in the large N 2 
limit |10] as 

Z M [N 2 , N,} ~ N^N^ exp(- 7M ^£) exp(/i c iV 2 + /if jV\), (3.3) 

v iV2 ' 

where the exponents are given for the pure gravity by au = —5/2 and 
Pm = 1/2- The non-universal constants jm, He and /xf depend on the 
choice of manifold conditions employed for the model according to the types 
of graphs or diagrams included. The matrix model partition function is then 
written as 




x exp{-(^-^ c )JV 2 -(^ B -^f)JVi}. (3.4) 

Here, we note the constraints in the sums over (N 2 ,Ni) are slightly relaxed 
owing to the strong damping factor exp(— ^mNI/N^. Sums in the partition 
function are approximated to integrals in the continuum limit by making 
scale factors for the area and the boundary length, A , and lo, tend to zero: 

Z[\,\ B ]^ [ X1 dx[ Vl dyx a yPexp(-—)exp(-\x-\ B y), (3.5) 

where we keep the regions of integrals finite in order to study behaviors 
toward the continuum limit, (xo,?/o) — ^ 0, as well as the large size limit, 
(£1,2/1) —7- 00. The cosmological constants are shifted and rescaled accord- 
ingly as A = (fj, - n c )/(jAo), and X B = {fi B - /if )/ (7/0). 
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We generalize the exponents (a, (5) through employing values obtained by 
the boundary Liouville theory, a = —1 — 1/b 2 and f3 = — 2 + 1/b 2 , in order 
to take into account matter fields, which are assumed to couple conformally 
with the space-time. We should remark that the difference in exponents 
Pm = 1/2 and f3 = —1/2 for the pure gravity is due to the way of counting 
distinct configurations. The matrix model respects the point group multi- 
plicity factor, Ni, in the state number counting. On the other hand the 
boundary Liouville theory regards those configurations which are members 
of the same point group symmetry are identical, since the universe stands 
by itself and does not have any reference frame. We shall employ the latter 
viewpoint in the followings. 

The integration over x converges, both in the limit xq — > due to the 
factor exp(— y 2 /x) for y > y > 0, and in the limit x\ — > oo because of 
exp(— Xx) for A > 0, to 

dx x a exp(-— ) exp(-Ax) = 2\ v/2 y- v K v [2y/\y], (3.6) 



Jo 



10 x X 

where we define v = 1/b 2 = — (1 + a), and where K v [z] is the modified Bessel 
function. We note that according to the small z behavior of the Bessel 
function, K v \z\ ~ z~ v , the integral converges even at A = 0, if we recall 
v > 1. On the other hand the integration over y, 

[ V1 dy y^exp(-\ B y)y-"K u [2V\y], 

requires a UV cut-off, y , while in large y regions it is possible to take the limit 
yi — > oo, if 2y/\ + X B > 0, since asymptotic form of the integrand contains 
the exponential factor exp{ — (2\/A + \ B )y}- The critical line separating the 
static phase and non-static phases locates along \ B + 2v^A = in the lower 
half-plane of the parameter space (A, \ B ), which we named as Bl, and along 
the line A = in the upper half plane as B3. 

Since all the statistical moments we are interested in are obtained by 
making use of the relation 



(—']) n + m f) n 
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where the statistical average is represented by brackets (...), we shall pay our 
attention on diverging parts of the integral which dominate in the continuum 
limit, yo 0. Utilizing the power series expansion of the modified Bessel 
function K v [2y/\y], and expanding the exponential factor exp(— \ B y), we 
obtain the diverging parts of the integral as 



\n+m 



C-lV 

I i • 1 1 V • 



n,m 



with 

2n+m— 1 — u 

T(v - n) 7-. 77, t \ ifz/>2n + m-l /n n . 

V > u+l-{2n+m) (3.9) 

if v < 2n + m — 1 , 



J (n,m) 



for v > 1, and for z/ = 1 we have 
1 A^ 



ZA.A' - 



I{(A B ) 2 - 2A + 4 lE X + 2A log A + 2A \ogy } logy ,(3-10) 

where 7^ represents the Euler gamma. 

From the partition function we can calculate statistical averages (x), (y), 
and correlation functions (A™ A™), of fluctuations, A x = x — (x) and A y = 
y — (y), for any pair of non-negative integers (n, m). For example, (x) and 
(y) are given by 

{x) = __ log z~-_ 

<H^, if-i (311) 

2(y logy ) if v = 1, 



('/) = ' T^TjlogZ 



<9A* 10gz Z {00 ) 



yo- (3.12) 



v 



They are independent of the parameters A and A , and the collapsing phase 
approaches smoothly toward the Bl and B3 boundaries. Since the average 
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values (A) = (x)/^ and (/) = (y)/j vanish in the continuum limit y — > 0, it 
is the reason we call this phase as collapsing. 
The second moment, (A 2 ,) is given by 

,2 - r~ ^2 



<A*> = 



d 2 



\ogZ 



v (20) 



£(10) 1 ' 
£(00) J 



<9A 2 *"°" Z (00) 

.4 " 



= < 



(u-ir 



~m lo S2/o 




(3.13) 



if i/ = 3 
if 1 < v < 3, 

where the dominant contribution of the last case will come from the non- 
diverging integral. In the same manner we can obtain other second order 
moments as, 



d 2 



d\d\ B 



fogZ 



"(oo) 



£(10) £(01) 
^(00) 



{-L.--^±i- }-hi y 3 ifz/>2 

-3^ log i/o ifi/ = 2 

if 1 < i/ < 2, 



(3.14) 



<9 2 



B\2 



log£ 



£(02) r£ 



'(00) 



(00) 



-2/o lo g2/o ifi/ = l. 

In the continuum limit all the correlation functions also vanish as 

,2 n +m\ if z/ > 2n + m - 1 

In + m — 1 
ifl<z/<2n + m — 1. 



(3.15) 



(A"A m 




y 2n+m logyo if* 



(3.16) 



4 Overview of non- stationary solutions of the 
master equation 

In order to investigate solutions of the master equation for non-stationary 
phases we shall analyze time evolution of statistical averages of moments 
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(x n y m ), instead of solving the partial differential equation. The master equa- 
tion provides equations of the motion as 

d\ X V ) „/„ i\/,n-2„,m\ i iN/„.n„,m-2\ 



dt 



= n(n - 1) (x n - 2 y m ) + m(m - 1) {x n y m ~ z ) 



where S = S(x,y) has been given by eq.(2.9). When the parameters (A, 
X B , j3) are chosen outside the region where the solution converges to the 
equilibrium distribution, it is beyond the case where standard procedures of 
the statistical mechanics are applicable. We shall require the master equation 
with the same function S(x, y) holds in the whole parameter space because 
of the equal probability axiom. 

Equations of the motion for (x) and (y) are given by 

^ = -A + (4> +a( i>, (4.2) 
dt x A x 

^ = -A B -2(^) + /3(i). (4.3) 
dt x y 

Since x and y correspond to the rescaled area and the boundary length, 
average values (y 2 /x 2 ) and (y/x) should be positive definite, even if (x) and 
(y) grow to infinity. Assuming (y 2 /x 2 ) — > c 2 and (y/x) — > Ci, eq.(4.3) can be 
written asymptotically as§ 

„ _ A + C2 + « ( I ); (4.4) 
dt x 

^ ~ -X B -2c 1 + P( 1 -). (4.5) 
dt y 

In these equations (l/x) and (l/y) are expected to behave like l/(x) and 
l/(y), respectively, when (x) and (y) get large. On the other hand, they 
will increase monotonically as (x) and (y) move toward zeros. Recalling the 



5 In the following 'tilde' (~) means asymptotically equal. 
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parameter a is negative, (x) increases indefinitely when —A + c<i is positive, 
while it goes to zero otherwise^ 

Therefore, when (x) — > 0, (y) should also tend to zero, which is the case 
we have already seen as the collapsing phase. On the other hand, if (x) — > oo, 
the increasing rate should be at most proportional to the Markov time. In 
fact, eq.(4.2) gives the asymptotic solution (x) ~ (—A + c 2 )t. 

The evolution of (y) depends on the parameters X B and (3. When /3 
is negative, d{y)/dt is negative and the solution becomes unstable except 
— — 1, which will be discussed in the later section as the polymer state. 
On the other hand when (3 is positive, (y) converges to an attractive fixed 
point for — A B — 2ci being negative, otherwise it expands with the maximum 
rate, as far as the parameter A lies outside the region of the collapsing phase. 

In summary, when X B < —1c\ and A < C2 , both (x) and (y) increase 
linearly in time, which corresponds to the expanding phase(P£ ). Asymptoti- 
cally c 2 and Ci tend to (y) 2 / (x) 2 and (y) / (x), respectively, and the critical line 
is located along vA = —2X B in the parameter space. The rest of the lower 
half-plane of the parameter space is occupied by the collapsing phase(.Pi ). 
In the upper half-plane of the parameter space (y) converges to a fixed point, 
and C2 and C\ go to zeros when A — > 0. This is the elongating phase(P<5' ), 
while the remaining region in the upper half-plane corresponds to the col- 
lapsing phase. 

These three phases are separated by three boundaries:the boundary Bl 
between PI and P2 is located along the line a/A = — 2\ B in the fourth 
quadrant. Since d(x)/dt is negative on this border, no physical solution 
exist, and the transition between two phases is abrupt, the boundary B2, 
which lies between P2 and P3, is specified by —A > and \ B = 0. Since 
d(x)/dt = —A and d(y)/dt = /3(l/y), there exists a solution where the space 
expands as (x) ~ 0(t) and (y) ~ (9(t 1//2 ), which we call the boundary state. 

4 Precisely speaking, because of the repulsive nature of the fixed point located at —A + 
C2 + a(l/x) = 0, the distribution bifurcates into the collapsing phase and the expanding 
phase with appropriate ratio, which tends 100% collapsing as —A + c-i approaches to 0. 
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The boundary B3 lies between P3 and PI, where A = and \ B > 0. Since 
d{x)/dt is also negative on this border, no physical solution exist, and the 
transition between two phases is sudden similar to Bl. 

In the following sections we shall show explicit forms of asymptotic solu- 
tions for these phases and boundaries. 



When we choose the parameters A and \ B negative, configurations with more 
triangles and more boundary edges become favorable, and the universe ex- 
pands with its maximum growth rates proportional to the Markov time. We 
have called this case as the expanding phase after the numerical simulation. 
In such a case it is common to separate fluctuations from their means as 
x — (x) + A x and y — (y) + A y . Equations of the motion are constructed for 
the averages (x), (y), and various moments of fluctuations (A™ A™), which 
are summarized in the Appendix B. Since the mean values, (x) and (y), in- 
crease with their maximum rates, we can expand equations of the motion in 
terms of the relative fluctuations, £ = A x / (x) and rj = A y / (y). As we shall 
see in a moment, the mean square fluctuations, (A 2 ,) and (A 2 ), also show lin- 
ear increase in time, the order parameters of expansions for x and y motions, 



y (O and y (i] 2 ), are proportional to t -1 / 2 . Therefore, solutions obtained by 

keeping the leading orders of expansion become exact asymptotically. 

For example, asymptotic equations for (x) and (y) are obtained from 
eq.(4.2) and (4.3) as 



5 The expanding phase (P2 ) 





d(x) 



d(y) 



dt 



dt 




(5.1) 



(5.2) 



Solutions of these equations are easily obtained by writing 




v x t, (y) ~ v y t, 
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where we keep only the leading terms, and require that v x and v y should 
satisfy 

v 2 

v x = -A + -|, (5.3) 
—X B - 2\ (5.4) 



v 



v 



The positivity of growth rates, v x and v y , defines two boundaries of the 
expansion phase in the 2- dimensional parameter space. The conditions, v x = 
and v y = 0, with keeping the ratio v y /v x {= r ) finite, give the boundary 
Bl along 4A = (\ B ) 2 with \ B < 0. This boundary coincides with that of 
the collapsing phase in the lower half-plane as we have seen in Sec. 3. There 
exists another boundary B2, where v x = —A and v y = 0, which lies along 
X B = with A < 0. 

Asymptotic equations for the 2-point correlation functions are written as 

d(*l) o ft) 2 f(*Z) (^) \ (5 . 5) 



dt (x) 2 ^ (x) (y) 

dJAaAy) (y) 2 r (A^) (A x A y ) 

dt ~ (X) 2 V (y) ( X ) 

They also grow linearly in t. Writing 

<A2)~ C (2,0K (A x A s >~c(l,l)t, (AJ> ~ c(0, 2)t. 
their coefficients are obtained to be 

2(4 + v x ) 



c(2,0) 

<± T 

c(l,l) = — -^J , (5.8) 

c(0,2) = 



4 + 4rg 
2(4r + v x ) 



4 + Ar 2 + v x 
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The asymptotic equation for the (n, m)-moment is separated into two 
branches, n + m = even and n + m = odd. For n + m = N with an even 
integer N there are N + 1 equations, 

d(A n A m ) 

K X dt y ~ ~ «(«-l)(Ar 2 A™> + m(m-l)(A^- 2 ) 

+ 2m— ( — s- T\)- 5 ' 9 

Inserting the asymptotic form of the solution, 

<A^A™> ~c(n,m)r^, 

into eq.(5.9), we obtain 

— — + °- + )c(n,m) 

v 2 v x v x ' 

= n(n — l)c(n — 2, m) + m(m — l)c(n, m — 2) 

+^ c ( n _l )m + l) + ^ c ( n + l ? m _i) } (5.io) 

fx fr 

which can be solved step by step starting with the n + m = 2 components 
with the initial input c(0, 0) = 1. 

For the n + m =od<i integer branch the equation has extra contributions 
from n + m =even integer components, 

n M!^ 3 (Ar 1 A™)-(Ar 1 A™)(A^) 



_ i (AgA-+ 1 )-(Ar 1 A-)(A a A,) 

( x ){y) 

(Ar 1 A^ 2 )-(Ar 1 A-)(Ag) , 

(y)/(Ar 2 A™- 1 )-(A«A- 1 )(A2) 



(x) V <*) 



2 



(A; +1 A-)-(AgA- 1 )(A x A,) A 

j ' (5 - n) 
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which should be added to eq.(5.9). Inserting the expected asymptotic form, 

into eq.(5.9) with eq.(5.11) we obtain 

/n + m — 1 2nr% 2m\ . . 

(—2— + ir + -) cKm) 

= n(n — l)c(n — 2, m) + m(m — l)c(n, m — 2) 

2wo 2mr 
H c(n — 1, to + 1) H c(n + 1, m — 1) 



+— [3r'{c{n+ 1,to) - c{n - l,m)c(2,0)} 

— 4r {c(n, to + 1) — c(n — 1, to)c(1, 1)} 

+c(n - 1, to + 2) - c(n - 1, to)c(0, 2)] 

— [r {c(n + 2, to — 1) — c(n, m — l)c(2, 0)} 

v x 

-c(n + 1, to) + c(n, to - l)c(l, 1)], (5.12) 

This equation can also be solved iteratively starting from n + to = 3, with 
initial inputs c(l,0) = c(0, 1) = 0, together with n + to = even solutions 
previously obtained. 

Near the boundary 51 of the expanding phase there are two cases of 
limiting procedures: i) take the t — > oo limit before taking (f^,^) — >■ 0, 
or 2) search a new state just on the boundary before taking the asymptotic 
limit. In the first case the expansion in terms of the relative fluctuations, £ 
and rj, is considered to be valid. As the parameters get close to the boundary 
values, the leading terms of expansion rates for (x) and (y) tend to vanish as 

A v y ~ -, , x 0, 



x 1+A ' y 1+A 

for a small 5 = — (\ B + 2\/A). On the other hand the growing rates for 
second moments stay finite towards Bl seen in eq.(5.8) for v x — > . Higher 
moments of correlation are obtained by solving eq.(5.10) and eq.(5.12), and 
then take the v x — > limit. 
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For the second limiting case we search the following type of solutions on 
the boundary: 

(x) ~ x ± t p + ..., (y) ~ Vl t p + ..., 

where < p < 1. The correlation functions we look for are considered to 
have the form, 



, v n-\-m I -| _ 

(A n A m ) ~ { c ^ m ) t ~^~ ~ if n + m = even 
1 Ci(n,m)t 2 i{ n + m = odd 

with < q < 1. At this stage we assume equations of the motion can be 
expanded in terms of £ and 77, and we take up to the second order, whose 
validity will be checked later. Inserting the asymptotic forms, equations for 
the averages (x) and (y) are given by 

p Xl t^ ~ «r*(i + ^V*) 

Xi v xf 7 

xf v xf xiyi yf / 

p Vl t^ ~ ^(l + ^V 2 ^ 

- 2^1|(#-%*}, 
xi v xf xiyi ' 

and for the 2-point correlation functions, 

gc 1 (2,0)^ 1 ~ 2-2a^%^^ 

x x 

Xi^ xf xiyi ' 

v xi 7 xiyi 

v xf xiyi 7 
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v? 

^ ^(0,2) CjiXV s 
xi ^ y\ xiyi ' 

There are three possible pairs of exponents (p, q), 1) (2/3, 1), 2) (1/2, 1), and 
3) (1/2, 1/2). The cases 1) and 3) give xi — y± — 0, and the case 2) does not 
have a solution due to negativeness of the parameter a. In summary, when 
the parameter 5 approaches to zero, the growing rates v x , v y get smaller, but 
nonetheless (x) and (y) increase indefinitely. On the other hand, as we have 
seen in the last section all the expectation values in PI tend to zero in the 
continuous limit y — )■ 0, and we conclude the transition at the boundary Bl 
is sudden. 

Let us turn our attention to the boundary B2, where the expansion rate 
for (y) vanishes as v y ~ (— A B )/(1 + 2v x ), while v x ~ —A stays finite with a 
negative A. In this case the ratio r = (y) / (x) approaches to 0, and coefficients 
of the two point functions tend to 

c(2,0) ~2, c(l,l) ~0, c(0,2) ~2- 



A + v x 

In general, coefficients of the correlation function (A"A™) behave for a small 
X B as 

c(even, even) ~ c(odd, even) ~ 0(1), 

and 

c(odd, odd) ~ c(even, odd) ~ (9(A B ). 

In the case when we set X B = before taking the asymptotic limit, we 
expect the following type of solutions, 

(x) ~ v x t + . . . , (y) ~ yit p + . . . , 

with < p < 1. Assuming that the expansion in terms of £ and 77 is possible, 
we have (x) ~ —At, while from the equation for (y), we obtain p = 1/2. 
Since the equation for (Ay) gives non-vanishing c(0, 2) = 2t> x /(4 + u x ), the 
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order parameter along the y-direction is \J {rf) ~ O(t ), which implies that 
the expansion is perturbative. We shall consider this case in the subsequent 
section. 

6 The elongating phase (P3 ) 

When A < 0, \ B > 0, and f3 > 0, the equation for (x), eq.(4.2), gives 
(x) — > v x t with v x = —A, assuming 



as (x) increases. Then, the equation for (y), eq.(4.3), has a stationary solution 
at the attractive fixed point, (y) = y±, satisfying 



where y\ is regarded to be a parameter at this moment. Since (A x ) also shows 
a linear increase in t, equations of the motion for the x-direction are expanded 
in terms of the relative fluctuation £, while equations for the y-direction are 
treated without separating the fluctuation. The set of equations of this type 
is named as the A-y scheme, and it is summarized in the Appendix C. 

In this section we choose (3 positive, otherwise d(y) / dt is always negative 
for \ B > 0, and the motion along the y-direction becomes unstable except 
f3 — — 1, which will be discussed in the next section. Keeping only the lowest 
order terms of expansion along the x-direction, equations of the motion in 
the asymptotic region are given by 




A B + /3(-) = 0, 



y 



d(A n x y m ) 



n{n - l)(A^- 2 y m ) + m(m - 1 + /3)(A> m - 2 > 



dt 



-— (A^-mA^A^- 1 ) 
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+ ^{(Ar 1 y m+2 )-(ArV")(y 2 )} 

+ ^- 3 {(K +1 V m ) - (ArV)(A 2 )}], (6.1) 

where the last three terms in the bracket [. . .] are reserved for the odd n case. 
We start solving the n = component, 



d (y m ) 1 , m,..m-2v o^(i/ m ) 



, , ~ m(m - 1 + /3) {y m - 2 } - - mA a (y™- 1 ) . (6.2) 

The solution is formally expressed as, 

(y m )(t) ~ mr 2m/v * f z 2m/Vx {-X B {y m - 1 ){z) + (m - 1 + /?)(y m - 2 )(z)}<^, 

(6.3) 

where we choose (x) = i; x i. Carrying out the integration explicitly, the first 
two solutions are expressed as, 

(y 2 ) ~ 2(l + ^-A B y 1 ) r ^ X + 2/2, (6.4) 

Ha: 

(y 3 > ~ 3 {(2 + /3) yi -A%} ' 



1 + ^ 

Vx 



+ A B (l + /j-A\) -xy + y 3 , (6.5) 

where y2 and y 3 are initial values of (y 2 ) and (y 3 ), respectively. Since the 
expectation value (y m ) should be positive, the solution (y 2 ) is acceptable only 
when the parameters satisfying 1 +/3 — \ B yi > 0. On the other hand, for the 
solution (y 3 ) the positivity condition requires 1 + j3 — \ B yi < 0. Combining 
these conditions, y\ is fixed to be (1 + (3)/\ B '. Continuing this procedure for 
larger m, we will find that the allowed asymptotic solution, (y m ) — > y m , is 
required to satisfy 

(m - 1 + ,%m-2 - A B y m _i = 0, (6.6) 
which gives the asymptotic value of (y m ) to be 

fm=II(-xB-j- ( 6 - 7 ) 

fc=l 
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For an even integer n the (n, 0)-component equation, 

^~«(»-l)(Ar 2 ), (6-8) 

gives the solution 

(A n x ) ~ d n t n/2 , 

n/2 

d n = 2™/ 2 n(2A;-l), (6.9) 

k=i 

for n > 2, and rf = 1. Since y ~ 0(1) and A s ~ Oit 1 / 2 ), inserting the 
asymptotic form, 

(A n x y m ) = c(n,m)t n / 2 , 
into eq.(6.1), the leading term follows 

^c(n, m)t n/2 - 1 ~ n(n-l)c(n-2,m)t n/2 - 1 +m(m-l + (3)c(n,m-2)t n/2 
2m 

_ — c ( n> m )t«/2-i _ mA B c(n, m - l)t n/2 , 

which requires c(n, m) to satisfy 

= ( m - 1 + /3)c(n, m - 2) - A B c(n, m - 1). (6.10) 
It gives the solution to be 

c(n,m) = n(^Mn,0). (6.11) 

From eq.'s (6.7) and (6.9) it is written as c(n,m) = d n y m , which indicates 
the lack of correlation between motions along x and y directions. 

We now proceed to solve the n = odd integer case. Since a correlation 
function of odd moment is not necessarily positive, we cannot utilize the 
positivity condition as the n-even case, and it is desirable to know the bound 
on exponents in order to select the acceptable solutions. Making use of the 
Schwartz inequality, 

(A>™> 2 < (K n )(y 2m ), (6-12) 
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together with the previous result for {(A^ n )} and {(y 2m )}, we have the 
bound, 

\(A»y m )\<0(t n ' 2 ), 

which is required asymptotically. 

Equations of the motion are given by eq.(6.1) with the additional terms 
reserved in [...]. Writing the asymptotic form of solution for the even n case 
as 

(A n x y m ) = d n y m t n /\ 
the contribution from additional terms is 

—d n+ iy m t 2 + —d n -i{y m+2 - y m y2)t 2 + — (d„+i - d„-id 2 )i/ m t 2 . 

(6.13) 

Since it has a different t dependence in the m = case from other cases, we 
shall solve the equation for m = and m^fl separately. 
Firstly, we set m = in eq.(6.1), and obtain 

^ ~ n(n - l)(Ar 2 ) + ^{( A " +1 ) " (AT'XA^)}. (6.14) 

This equation can be solved iteratively starting from n = 3 using the solutions 
for even n moments, resulting 

(A") ~ dJ n -V/ 2 \ogt, 
~ 2n{n-\) ~ 

d n = —, -—d n -2, (6. 15) 

(ra-3) 

for n > 5, with d^ = (24a /v x ). Next, we consider the m^O cases, and we 
get the same relation as the previous even n case: 

(m - 1 + f3)(A n x y m ' 2 ) - A^A^T" 1 ) ~ 0. (6.16) 

Solutions of this equation were already known to be 



(A*y m ) ~ (A;>n 



k=l 

rfW {n - 3)/2 logt, (6.17) 
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which also shows the lack of correlations between x and y motions. 

Approaching toward the boundary B2 from above, \ B — > + , the mo- 
ments, (A"?/ m ) diverges as 0((\ B )~ m ). In order to understand the nature of 
transition between the expanding phase and the elongating phase, we need 
to know properties of the boundary state along B2, which will be discussed 
in the subsequent section. On the other hand (A™?/ m ) does not depend on A 
asymptotically, and the transition at the boundary B3 is abrupt according 
to the limit in the neighboring phase PI. 

7 The polymer state and the boundary state 

There exist two special states for certain values of the parameter set (A, \ B , 
{$). When (3 is zero or negative, d(y)/dt is always negative for X B > 0, and 
(y) tends to be unphysical. In the negative /3 region there exists a special 
case locating at j3 — — 1, where not only y 2 but also all {y m } are according 
to eq.(6.7). In this case the distribution function p(x,y;t) approaches to 



We may call this case as the polymer state. In this state only the motion 
along the x-direction evolves as (x) — > v x t with v x = —A, and 



except the n = 1 component, i.e. (A x ) = 0. The coefficients {d n } and {d n } 
are the same as those of the elongating phase given by eq.(6.9) and eq.(6.16), 
respectively. 

There is another special state which appears along the boundary B2. For 
an even integer n the equation of motion, 



p(x,y;t) ->p(x;t)6(y). 




(7.1) 



d{A%y m ) 



n(n - l)(A™- 2 y m ) + m(m - 1 + /3)(A> m ~ 2 ) 



dt 



2m 



(Ky m ), 



(7.2) 



(x) 
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has the solution with the product form, 



where each component obeys, 

n(n-l)(AT 2 )- (7.3) 



^ ~ m(m-l + /3)( 2/ - 2 )-^( 2/ m ). (7.4) 

The first equation is the same as the elongating phase and we already have 
the solution (A£) ~ d n t n / 2 with the same d n given by eq.(6.8). For the second 
equation writing (y m ) ~ y m t m ^ 2 , we obtain 

4 _ 

y m = 2(m - 1 + /9)j/ m _ 2 (7-5) 

For an even integer m we solve it iteratively with the initial value yo — 1, 
and we get 

2 m m / 2 

^ = (tta) 2 n ( 2A; - 1 + ( m ; even )- ( 7 - 6 ) 

~r Ua! fc=i 

In order to obtain solutions for an odd integer m we need to give a value for 
yi, which is regarded as a parameter at this moment. Then, the solution is 
given by 

O m^l ("»-l)/2 

Vm = (j— x) 2 II ( 2A; + (m : odd). (7.7) 

Dl fc=l 

For an odd integer n we need to add extra-terms coming from the even 
n solutions in the equation of motion, eq.(6.1). Since contributions of these 
extra terms for m = and m ^ have different properties similar to the 
previous case, we treat them separately. For m = the equation is same as 
the elongating phase and the solution is known to be 

(A n x ) = dJ n -W 2 \ogt, (7.8) 

with the same coefficient d n as eq.(6.16). 
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For m 7^ case the additional terms are written by using the n-even 
solutions, (A™y m ) ~ d n y m exp{ — (n + m)/2}, to be 

2777- 

vl v 2 x 

where we keep terms having the order ^( n + m - 3 )/ 2 . We note the additional 
contributions can be rewritten as 



3 + ^ 



m 

dn+iyini 



using the relationships for d n , eq.(6.9). Writing the asymptotic solution as 

except the (1, 0)-component, which is by definition, the equation of motion 
is expressed as 

ra + ra — 1 - , lN - „ , , m ? - 
2 «nl/m = n(n - l)d n _ 2 y m + m(m - 1 + fj)d n y m _ 2 

- Qnl/m + — -f -ofln+iym, (7-9) 

w - 1 + ^ u ' 

except the equation for the (n, m) — (1, 2) moment, which involves the (1, 0)- 
component, 

4 

^ 1 + ^ u z 

Imposing 

d n = 2nrf„_ 2 , (n : odd) (7.11) 
with di = 1, eq.(7.9) becomes 

4 3 + -2- 4 _ 

y m = 2(m - 1 + /3)y m _ 2 2/ m + : ■ "4 -zy m , (7-12) 

^ 1 + ~ x V x 

where we use the relation 

d n +i/d n = d n -i/d n -2 — d2/di = 2. 



dm = dm + v -f^d 2 y2- (7.10) 
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We note eq.(7.10) becomes the same as the m = 2 component of eq.(7.9) 
when we set y = 0. 

The solution of eq.(7.11) is given by 

(n-l)/2 

d n = 2 {n ~ 1 ^ 2 J] (2k + 1) (n: odd), (7.13) 
fc=i 

while eq.(7.12) can be solved for even m and odd m separately, using the 
solution y m given by eq.(7.6) for the even m case and eq.(7.7) for the odd m 
case. For an even m we solve it iteratively with the condition y = 0, while 
for an odd m the solution has the same form as the even case with the initial 
input value y\, which we assume to be a parameter at this moment. We give 
here the solutions without giving details of derivations, 



>(t^) 2 T\T=i(2k- 1 + P) if m = even 

V - 7 r^l (7.14) 

:^ayi + yi)(iqV) 2 n& 1)/2 (2A; + /3)yi if m = odd, 

where we write a = 4(3 + ^)/v 2 x /{l + ^) 2 . 

Until this point we obtain all quantities of our interest in the asymptotic 
limit with two initial input parameters, y~\ and y\, which correspond to (y) 
and (A x y), respectively. Since their equations of the motion, 

f> ~ 2 M (, 15) 
<JM. . , (Ai ,- V2 <^ + W_|M, ( , 16) 

involve moments with the negative exponent, it seems we are eventually 
forced to solve infinite towers of equations with negative integer arguments. 

In order to avoid these difficulties we shall make use of the following 
relation for (y -1 ), 

i oo n 
\y/ n=0m=0 
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Vi 



n=0 m=0 




i oo n 

TEH)1 " ]gr 1/2 - (7.i7) 



Since y m is written as eq.'s (7.6) and (7.7), eq.(7.14) indicates (y x ) is ex- 
pressed by £/i. From this relation and the eq.(7.12) we have the equation, 

I(i + 1)^ = A v v (-ir i n 1 — , (7.i8) 

which can, in principle, determine the parameter y~\. Similar analysis can be 
applied for (A^y -1 ), giving 

z ^ Z/i „ =0m =o \m J Vi v x 

Again, y m is known from eq.(7.11) up to two parameters, y~\ and yi, and the 
remained parameter y\ should be unambiguously determined. 

This concludes analyses of solutions of the master equation obtained 
through equations of the motion for all the possible moments, and they 
show the asymptotic solutions are exact and are determined without any 
ambiguous parameters other than (A, A B , pi). 



8 Discussion: Evolution of the universe 

After studying statistical properties of the 2-dimensional triangulated man- 
ifold, we are ready to introduce the space and the time as the frame of 
reference. It is natural to choose space and time coordinates according to 
the global property of the topology. The manifold we have constructed has 
the two-dimensional disk topology. Peripheral boundary of the disk fits to 
the space coordinate from the symmetry of forward and backward direction. 
On the other hand the direction perpendicular to the boundary is not sym- 
metric, and we can easily recognize the inside as past and the outside as 
future when it is regarded as time. 
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Among the solutions of the equation of motion the average area, A(t), 
and the average boundary length, Z(r), are the two fundamental quantities 
for the space and time. Their evolution are described by the length of Markov 
chain, which we write r, here. We consider these two quantities are related 
each other through the 'time' t, which we may call the physical time, as 



where c takes into account of adjusting scales of the physical space and the 
physical time. 

The time we defined in eq.(8.1) is a function of r, which is determined by 



From the solutions, (x), and (y), we have the area and boundary length 
A = (x)/7, and I = (y)/7, and we can calculate the relationship between t 
and r, for each phases as the followings. 

PI: A(t) ~ Z/0/7, and /(r) ~ yo/7 for v > 1 




(8.1) 




(8.2) 



ct = (no running time). 



We note that the crunch occurs as a result of collapsing. 
P2: A(t) = v x t/^, and Z(t) = v y r/^ 




P3: A(t) = u x r/7, and Z(r) = yi/7 



ct = — t, and l(t) = — (stationary universe). 
Boundary State: A(t) = v x t/^, and /(r) = yir 1 ^ 2 /^ 



ct = —\fr, and lit) 
Vi 




t ('Big Bang' universe). 
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We employ the name Big Bang as the power law expansion in comparison to 
the inflation as the exponential expansion. 

Polymer State: A(t) = v x t/^, and Z(t) = yi/7 

v 

ct = — r, and l(t) = (stationary universe). 
The polymer universe appears when /3 = — 1, where 7 _1 = 

9 Summary 

We have started creating a two-dimensional universe under the working hy- 
pothesis, There were no laws but rules at the beginning of universe. The 
rules we have proposed consist with the three main axioms: 1) existence of 
the elementary unit of universe, 2) existence of neighbors of a world, and 3) 
equal probability of appearance of worlds. Based on these axioms we have 
obtained the master equation which describes the emergence and evolution of 
a two-dimensional space-time. The equation reveals three prominent phases, 
i.e. the expanding, the elongating and the collapsing phase, and two signifi- 
cant states, named as the polymer and the boundary state. For each phase 
and state statistical averages of the area, the boundary length, and corre- 
lations between those variables and their fluctuations are obtained. Since 
these solutions are expanded in terms of the order parameter, which is the 
inverse square root of the Markov time, they are considered to be asymp- 
totically exact. Moreover, solutions contain no ambiguous parameters in the 
asymptotic limit except three basic inputs, i.e. the cosmological constant A, 
the boundary cosmological constant \ B , and the matter content j3. 

The physical time is introduced from the relation between the area and 
the boundary length, after assigning the boundary as the physical space. 
In terms of the physical time solutions of each phase show the characteris- 
tic behavior of space evolution, 1) inflationary expansion in the expanding 
phase, 2) stationary universe in the elongating phase, 3) point universe in 
the collapsing phase, 4) stringy universe in the polymer state, and 5) linear 
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expansion in the boundary state, which we have called the Big Bang expan- 
sion. The variety of solutions indicates the ability of describing most of the 
essential features of the universe even in the 2-dimensional space-time. We 
remark all the possible type of time appeared in the two-dimension seems to 
be universal for any manifold with one boundary according to the numerical 
simulation for the 4-dimensional case. 

In order to investigate the space-time evolution in higher dimensions we 
need to achieve method to enumerate triangulations. Because of the lack of 
analytical means numerical efforts to obtain some hints for functional forms 
of the distribution function are unavoidable. The dynamical triangulation 
with boundaries has been performed sometime ago for searching continuous 
transition [161 E] in space-time dimension 3 and 4. Although our purpose 
is to enumerate triangulated graphs of open manifolds, those investigations 
suggest possibikities to apply the dynamical triangulation prepared for closed 
spaces even for open cases. 

To accomplish the creation of universe it is the most important to answer 
following two questions among many. Firstly, how the three parameters 
are determined at the very beginning? Part of this question will be solved 
by moving up to higher dimensions. In a four- dimensional universe, where 
the conformal symmetry is broken and the cosmological constants are not 
constant any more, the effective cosmological constants will be allowed to 
make time evolution. In the 2-dimension the cosmological constants obtained 
by the matrix model are large, and the universe is in the expanding phase. 
If the conformal symmetry is broken, these two parameters are not constant 
anymore, and they will converge to the small fixed point in the parameter 
space. For the matter part of parameters we may say that any combination 
of fields is possible. In case we are asked to explain the raison df?tre of 
the matter content of this universe, we need to accept the anthropological 
principle at this moment. 

Second problem is the birth of quantum dynamics, which we think the 
primordial dynamics of universe. In contrast to the master equation which is 
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the same type as the diffusion equation, dynamical equation is non-diffusive. 
It is the hardest part of problem to answer how the complex function becomes 
inevitable in the course of evolution of space and time. Unfortunately, we do 
not have any answer to this question. 

In spite of the lack of quantum dynamics we consider the creation of space- 
time can stand by itself. Same kind of situation arises in evolution of metals: 
firstly a metal lattice is formed which is dominantly the classical statistical 
process. Later, conduction of electrons and the electron-phonon interaction 
emerges in the metal which are described by quantum mechanics. Afterall, if 
the evolution of universe at its beginning is inflational, the quantum dynam- 
ical correlation cannot influence the space-time dynamics beyond its limited 
event horizon, and the evolution of space becomes inevitably statistical. 
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Appendix 

A The master equation in the asymptotic limit 

The Markov process in the simplicial space is described in terms of the de- 
screte variables (N 2 ,Ni) and the descrete time n: 

p a [n + 1] - p a [n] = J2{Pb[n\w b ^a ~ PaH^i.}, (A- 1 ) 

6(a) 

where w a ^ = {1 + exp(Sb — S'a)} -1 with 

/V 2 

S a = S(N 2 , N,) = -a log N 2 - f3 log N, + 7 -f + /iiV 2 + ^N,. (A.2) 
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Here, the sum over b(a) is restricted within the four neighbor points of a = 
(N 2 ,Ni) connected by the four moves, Mi : (N 2 + 1,N ± + 1), M 2 : (N 2 - 
l,Aq + l), M 3 : (iV 2 + l,iVi-l), M A : (N 2 - 1, Ni - 1). Introducing the three 
scale parameters (A ,l ,r ) for the physical variables, the area A = A N 2 , 
the boundary length / = IqNi, and the Markov time r = r n, the probability 
distribution function is translated as 

P(N 2 ,N 1 )[ n \ ~p(a,1;t). 

In the continuum limit, where (A , l , r ) — > 0, the finite shift of a veriable of 
a function, like f(x + xo) = exp(xod/dx)f(x), is expanded as {1 + x^d/dx + 

(l/2)(x d/dx) 2 H }f(x). In the l.h.s of eq.(A.l) the lowest order surviving 

term is C(r ), while the r.h.s of eq.(A.l) the expansion starts with the second 
order in A and / . We shall keep only those leading terms for the continuum 
limit assuming O(r ) ~ O(A^) ~ 0(Iq) to obtain 

d r 2 , d 2 dS a d d 2 S a , { d 2 dS a d d 2 S a , x 

ro^ P [A,l;r] = { Afi (—+—— + )+/„( + + )j p [A,Z;r]. 

(A.3) 

It is free to choose the scale parameters so that they satisfy r = A^ = 1%. 



B The A-A scheme 

When the two parameters A and X B are chosen in the region of the expanding 
phase, which locates mostly in the third quadrant and partially in the fourth 
quadrant, the area and the boundary length increase at their maximum rates 
in time. Variables x and y are then separated to mean values (x), (y) and 
fluctuations A x = x — (x) and A y = y — (y). Since the second moments of 
fluctuations (A^) and (A^) also increase proportional to time, the relative 
fluctuations £ = A x /{x) and i] = A y / (y) are of the order Q{t~ l l 2 ). The 
equation of motion for averages of x, y are written as 

d(x) (i/) 2 , (l+>)) 2 , , a 1 

- - A+ wW + wW' (B ' 1) 
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and moments of fluctuations (A"A™) as 

[ ' dt '"' = »(» - 1 ) <AT 2 A™} + m(m - 1 ) (AJAJ'" 2 ) 

+ n {x yA a * A »t( 1+? )2 <(l+?) 2 '- f 

- 2ro |>(A;A r {i±f-(i±|>}). (B.3) 

C The A-y scheme 

When the boundary cosmological constant changes its sign (i.e. X B > ), the 
boundary length stops increasing and the universe elongates as time goes by. 
In this case expanding the equation for fluctuations along the y— direction 
cannot be justified, while the motion along x-direction keeps expanding with 
its maximum rate. Then, we shall consider equations of the motion without 
splitting < y > and < A y , and we obtain the equation for the average (x), 

dt - A+ (x}* { (i+o> )+ (c - 1} 

and the correlation (A"|/ m ), 

d(Kv m ) 



dt 



= n(n - I) (A™~ 2 y m ) + m(m - 1 + /3) (A> m ~ 2 ) 



2m^(A; s "' r ^> - ™A»(Ay-'). (C.2) 
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